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ABSTRACT 

We investigate a class of pulsar wind nebulae that show synchrotron emission 
from a thick toroidal structure. The best studied such object is the small radio 
and X-ray nebula around the Vela pulsar, which can be interpreted as the result of 
interaction of a mildly supersonic inward flow with the recent pulsar wind. Such 
a flow near the center of a supernova remnant can be produced in a transient 
phase when the reverse shock reaches the center of the remnant. Other nebulae 
with a thick toroidal structure are G106.6+2.9 and G76. 9+1.0. Their structure 
contrasts with young pulsar nebulae like the Crab Nebula and 3C 38, which 
show a more chaotic, filamentary structure in the synchrotron emission. In both 
situations, a torus-jet structure is present where the pulsar wind passes through 
a termination shock, indicating the flow is initially toroidal. We suggest that the 
difference is due to the Rayleigh- Taylor instability that operates when the outer 
boundary of the nebula is accelerating into freely expanding supernova ejecta. 
The instability gives rise to mixing in the Crab and related objects, but is not 
present in the nebulae with thick toroidal regions. 

Subject headings: ISM: supernova remnants — pulsars: general 



1. INTRODUCTION 

The structure of a pulsar wind nebula is determined by the interaction of the relativistic 
pulsar wind with its surroundings, which are initially the parent supernova and later the 
interstellar medium. Pulsar nebulae are best observed by thei r synchrotron radio and X-ray 



emission, and three classes of ne bulae can be characterized ( IGaensler fc Sland 120061 ) . The 



first, typified by the Crab Nebula ( Hester! 2008 ). shows a small toroidal structure close to the 
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pulsar, and a more symmetric, larger scale structure. The outer boundary is formed where 
the shocked pulsar wind expands into the freely expanding gas of the parent supernova. The 
second, typified by Vela X, is formed after the reverse sho ck wave from the int eraction with 
the interstellar medium moves back to the pulsar nebula (IBlondin et al.ll200ll ). This action 
displaces the pulsar nebula from th e pulsar, leaving a 'relic' nebula. The third, typified 
by the Mouse (IGaensler et al.l I2004J ) . shows a trail of emission as the pulsar moves rapidly 
through the interstellar medium. 

Here, we consider members of another class of pulsar nebula that is characterized by a 
thick toroidal emission region. In Section 2, we treat various members of this class, with spe- 
cial attention to the recently formed nebula surrounding the Vela pulsar. The distinguishing 
attributes of this class of pulsar nebula are discussed in Section 3. 



PULSAR NEBULAE WITH THICK TOROIDAL EMISSION REGIONS 



2.1. The Vela Nebula 



The best studied example of a remnant with a likely asymmetric reverse shock is the 
Vela remnant, where a strong case can be made that a reverse shock from th e N side of the 



remn ant has displaced the original pulsar nebula (Vela X) from the pulsar ( IBlondin et al. 



200ll ). While Vela X is now a relic nebula, the Vela pulsar has created a new nebula in its 
vicinity. Deep radio ima ges of the region aro und the Vela pulsar show evidence for a bow 
wave (Figs. 9 and 10 of iDodson et al.l l2003al . Fig. 1). The orientation of the bow wave is 
consistent with a flow from the N direction; the fact that the emission is strongest along 
a NE-SW line can be attributed to the orientation of the pulsar equatorial pl ane with this 



direc t ion. Deep X-ray images sho w similar struc t ure to that seen in the radio (jPavlov et al. 



20031 : Kargaltsev & Pavlov! 12004 ). IPavlov et all (120031 ) note that the X-ray jet to the NW 



tends to flop to the S, indicating the presence of a flow in that direction in the shocked pulsar 
wind. 

The presence of a weak bow shock in Vel a might seem surprisi n g bec ause the pulsar 
is centrally located in the larger remnant and Ivan der Swaluw et al.l (120041 ) showed that a 
moving pulsar must be a t least 68% of the way to th e outer shock in order to be moving 
supersonically. However, Ivan der Swaluw et al.l (120041 ) assume that the supernova remnant 
flow is described by the Sedov solution. There is a transition period before the Sedov solution 
is closely approximated. Considering a centrally located pulsar, the postshock region of the 
reverse shock (speed v s ) is characterized by a velocity 0.75u s and sound speed c s = (5p/3pY^ 2 , 
with p = 0.75po^s and p = 4p from the shock conditions, where po is the preshock density. 
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The result is a Mach number A4 = v/c s = 3/5 1 / 2 = 1.3 flow. If the pulsar has a velocity 
that carries it away from the center, it is comoving with the freely expanding gas, so that 
the properties of the postshock gas relative to the pulsar are the same as just described. If 
the reverse shock front has recently passed over the Vela pulsar, a mildly supersonic flow is 
naturally produced. 

The observed radio emission marks the shocked pulsar wind and thus the black line 
in Fig. 1 delineates the contact discontinuity (CD) between the shocked pulsar wind and 
the shocked surrounding medium. The half-opening angle of the CD at large distances 
is ~ 30° . The opening an gle of the bow shock is expected to be larger (e.g., Fig. 8 of 
van der Swaluw et al.l l2004j ). In flow past a blunt object, the half-opening angle, (3, of the 
shock wave that is generated obeys sin/3 = c s /v = where c q is the gas sound speed, v 



is the gas velocity, and M. is the Mach number (e.g.. lLandau &: Lifshitzlll959l ). The expected 
M. = 1.3 flow thus gives a half-opening angle for the bow shock of 50°, which is roughly 
consistent with the observations. The implication is that there is a mildly supersonic flow 
past the pulsar nebula. 



The direction of pulsar motion is known (position angle 301°. lDodson et al.ll2003bl ) and 
does not line up with the displacement of Vela X or the bow shock discussed here. The flow 
velocity involved in the displacement of Vela X is v a = 500(rd/5 pc)(td/2500 yr)" 1 km s _1 , 
where is the displacement distance estimated for the center of the relic pulsar nebula and 
i<2 is the a ge since the displacemen t began. The temperature in the central region of Vela is 
~ 1.2 keV (ILu fc Aschenbachj|2000l ). yielding a sound speed of 560 km s -1 , although there are 
also cooler components. The flow velocities associated with the reverse shock are apparentl y 
much larger than the projected velocity of the Vela pulsar, 61 km s _1 (IDodson et al.ll2003bl ). 
so that the pulsar velocity does not play a significant role in the orientation of the flow past 
the pulsar, in agreement with observations. If the surroundings of the pulsar have been 
cleared so that the pulsar interacts directly with the thermal gas, a mildly supersonic flow 
past the pulsar is expected. 

Some properties of the Vela pulsar nebula are given in Table 1. The nebular pressure 
P n at the termination shock (assumed to be the torus radius Rtorus) is given by pressure 
balance P n = E / (Att R^ orus c) = 2 x 10~ 9 dyn cm -2 for the parameters in Table 1, where E is 
the spindown power of the pulsar. If half of this pressure is supplied by the magnetic field, 
we have B = 1.6 x 10 -4 G. Assuming 7 = 4/3 for the shocked pulsar wind and a uniform 
pressure, the energy (within the nebular radius R ra dio) is E neb = 3.6 x 10 45 ergs. The 
minimum energy, E min , required to produce the synchrotron emission from the Vela nebula 
can be estimated. The total flux of the radio nebula is 1 Jy and has a flat spectrum between 
1.4 and 8.5 GHz (IDodson et al.ll2003al ). Extrapolating the low energy X-ray emission, with 
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ol = —0.66 where R, oc v a (IMangano et al.ll2005l ). yields a break frequency of 3 x 10 10 Hz. 
Equation (46) of IChevalierl ( 120051 ) for the minimum energy then yields E min = 4 x 10 44 
ergs. The energy estimated above is ~ 10 times the minimum. The time for the pulsar 
spindown to produce this energy is E ne b/E = 17 yr. Radiative losses are small and adiabatic 
expansion losses are not expected to change this estimate by a large factor (< 2), so the age 
of the particles in the radio nebula is relatively small. Older particles a re swept downstream 
in the SE direction, as can be seen in deep radio and X-ray images (IDodson et al.l l2003al : 
Kargaltsev fc Pavlov! UOM ) . 



As mentioned above, the spectrum of the to tal radio flux is flat, o r slowly declining, 
which is typical of pulsar wind nebulae. However, IDodson et al.l ()2003aJ) find large (several 
0.1 's) variation in a over the nebula when examined over the separate frequency ranges. 
Also surprising is that in some places a > 1 /3, which is not expected for optically thin 
synchrotron radiation. iBietenholz et al.l (119911 ) measured a = —0.6 between 1.4 and 5 GHz 
for the bright part of the radio nebula to the NE, which also indicates variations in spectral 
index. This situation contrasts with that of the Crab a nd related nebula e , whe re spatial 
variations in radio spectral index are small. For the Crab, IBietenholz et al.l (Il997| ) find that 
the a is constant to an accuracy of 0.01. Although the limit on variation is weaker in the case 
of 3C 58, the spectral index is also fairly constant ( IBietenholz et al.ll200ll ). These properties 
suggest that there are mixing or diffusive processes going on in the Crab and 3C 58 that are 
not present in the Vela nebula. 

The morphology of the Vela radio nebula contrasts with the pre-reverse shock pulsar neb- 
ulae. The magnetic fi eld is predominantly toroidal, with the linear polarization reaching 60% 
( IDodson et al.ll2003al ). Filamentary structure with a scale comparable to the termination 
shock radius, as seen in younger nebulae, is not apparent. The radio nebula is smaller than 
in the pre-reverse shock cases, but R ra diol ' Rtorus = 5 — 6; the angular size of the structures is 
comparable to the case of the Crab Nebula. We suggest that the absence of structure is due 
to a combination of the facts that the outer boundary is not subject to the Rayleigh- Taylor 
instability in this case and that only relatively recent effects of the pulsar wind are observed. 
The strong emission in the pulsar equatorial plane can be attributed to the stro nger pulsar 
powe r in this direction, as expected in a split monopole model for the pulsar wind (jBogovalov 
19991 ). The dynamics of the nebula depend on a complex 3-dimensional situation because 
the rotation axis is not aligned with the direction of the surrounding flow. The fact that a 
jet-torus structure is present around the puls ar suggests that the magnetization parameter 
a « 0.1 in the pulsar wind (jBucciantinill201ll ) and that the magnetic field reaches equipar- 
tition strength in the downstream flow. Once the magnetic field reaches this strength in a 
bow shock situation, it tends to increase the size of the wind nebula between the termination 



shock and the contact discontinuity (jBucciantini et al.ll2005al ). so the contact discontinuity 
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is at several times the termination shock radius, as observed in the Vela nebula. 



Another indication of significant magnetic forces in the nebula i s the central position 
of the pulsar in the inner X-ray nebula (Fig. 1, INg fc Romanil 120041 ) . In a hydrodynamic 
model in which the pressure is constant between the contact discontinuity and the pulsar 
wind termination shock, the ratio of the termination shock radiu s on the upstream side to 
that on the downstream side is (Ivan der Swaluw et al.ll2003l ). Assuming that the inner 
X-ray nebula represents the termination shock and Ai ~ 1.3, the displacement expected in 
the hydrodynamic model is not present. However, significant pressure forces from a toroidal 
magnetic field give a decreasing pressure with radius and the press ure gradient can be larger 



than the pressure differential from the bow shock estimated here ( IKennel fc Coronitilll984 
Emmering &: Chevalierlll987l ). 



We have examined simple models for the projected radio emission in which there is 
uniform emission from a toroidal region with an inner and out er radius, and an op ening 
angle. Fig. 2 shows a model that resembles the radio structure (IDodson et al.l l2003al ) : the 
outer radius is 5 times the inner radius, which is assumed to be at the position of the inner 
X-ray torus, the half-angle of the disk is 35°, and the angle between the line of sight and the 
axis of the torus is ( = 50°. The emissivity is assumed to be uniform and isotropic, which 
would be appropriate for a region with a tangled magnetic field. However, the observed 
radio polarization at 5 GHz ( IDodson et al.ll2003al ) is the maximum expected for a flat flux 
spectrum, suggesting that the magnetic field is ordered. Emission in a purely toroidal field 
would enhance the emission along the minor axis compared to the major axis and thus reduce 
the agreement with the observed morphology. Models for the synchrotron emission from the 
central region of the Cr ab based on relativistic MHD simulations show r easonable agreement 
with the observations (jBucciantini et al.ll2005bl ; iDel Zanna et al.ll2006l ). The outer parts of 
ring-like str uctures are enhanced w hen the magnetic field i s taken to be isot r opic a s opposed 
to toroidal (IDel Zanna et al.l 12006k their Figs. 3 and 4). IDel Zanna et al.l (120061 ) also find 
the effects of the mildly relativistic flow are significant. In our case, the radial velocity drops 
beyond the termination shock so that relativistic effects are not important. 

The observed em ission from the Vela nebula shows a high surface brightness to the NE 
(IDodson et al.ll2003al ). presumably due to the low Mach number flow. The flow is probably 
also responsible for the higher surface brightness to the NW along the minor axis, and pos- 
sibly for the brighter X-ray "outer jet" to the NW even though its motion is e xpected to be 



direc t ed away from us bas ed on the brightness distributions of the inner tori (jPavlov et al. 



20031: INg fc Romanil 120041) . The angle ( in our model is somewhat less than that found by 
Ng fc Romanil (120041 ) for the inner toroidal X-ray structure around the pulsar (( = 63.6°); 
again, the difference may be related to the bow shock flow. Three-dimensional MHD simu- 



-6 - 



lations will ultimately be needed to model the emission. 



The X-ray emission is observed to extend out to the same place as the radio (jKargaltsev fc Pavlov 



2004 ). The age of the nebula and the magneti c field strength indicate a synchrotron cooling 



break in the X-ray range and iMangano et al.l (120051 ) find an apparent cooling break in the 
spectrum at 12.5 keV. This would explain the similar extent at radio and X-ray wavelengths. 



2.2. G106.6+2.9 



Another nebula with similar properties to the Vela nebula is G106.6+2.9 (also known 
as the Boomerang), associated with J2229+6114 (Table 1). Unfortunately, the distance is 
poorly known, w ith estimates of 80 pc based on low radio depolariz ation and interactio n 
with neutral gas ( iKothes et al.ll200ll ). 3 kpc based on X-ray absorpti on ( jHalpern et al.ll200ll ). 
and 7.5 kpc based on pulsar dispersion measure (lAbdo et al.l 120091 ) . We scale the distance, 
d , with d% = dj (3 k pc). The Boomerang has considerable radio emission on one side and 



Kothes et al.l ( 120061 ) identified this structure as the result of an asymmetric reverse shock, 
although a larger scale thermal remnant is not observed around the nebula. In this case, the 
pulsar E and Rtorus suggest P n = 3.3 x 10~ 10 dJ 2 dyn cm~ 2 . The magnetic field and the radio 
nebula energy, computed as for Vela, are B = 6.5 x lO -5 ^ 1 G and E neb = 2.7 x 10 47 <i3 ergs. 
The minimum energy in the radio nebula is 2.4 x 10 46 g?3 7 ^ 7 ergs. The time for the pulsar to 
produce the nebular energy is again relatively short, ~ 390g?3 yr. 



The radio nebu la is highly polarized, up to 50%, indicating a toroidal field component 
( iKothes et al.l 120061 ). The radio spectrum sho w s a = —0.11 ± 0.05 below 4.3 GHz and 
a = —0.59 ± 0.09 above, which iKothes et al.l ( 120061 ) interpret as a synchrotron cooling 
break; they note that there is a steeper spec t rum f arther from the pulsar. To obtain such 
a low cooling break frequency, IKothes et al.l (120061 ) need a magnetic field of 2.6 mG and 
particle lifetime of 3900 yr in a model with d% = 0.27. We find that this field would produce 
too high a pressure in the nebula and the age is longer than expected, so that the spectral 
break is likely to be intrinsic to the injection. 



Kothes et al.l (120061 ) propose that the variations in rotation measure across the remnant 



are intrinsic to the remnant because of the apparent symmetry about an axis defined by the 
polarized emission. For the magnetic field we have derived, an electron density n e > 1/f 
cm -3 , where / is the filling factor would be needed in the remnant to produce the rotation 
measure. In our view, n e is much lower than this value because the pulsar wind bubble has 
smoothly expanded into the surrounding medium, so the apparent alignment of the rotation 
measure pattern is a coincidence. 
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2.3. G76.9+1.0 

The remnant G76. 9+1.0 has a similar radio morphology to the Vela nebula, and Chandra 
observations show t hat the radio lobes have a similar orientation compared to the axis of 
the pulsar rotation (jArzoumanian et al. 2011 ). However, it differs from the Vela nebula in 



2 important ways: it is much larger (29 x 35 pc at a distance of 10 kpc) and the large scale 
nebula has symmetry about the pulsar. These properties indicate that the pulsar nebula 
has not been affected by the reverse shock wave. We suggest that the inward mixing of 
thermal filaments has not been so efficient in this case so that a toroidal emission region has 
a chance to build up, as in the case of the Vela nebula. The instabilities are expected to be 
strongest during the constant E phase of the pulsar evolution. Once a pulsar is no longer 
adding substantially to the energy of the remnant, the acceleration ceases. This might be 
expected in a large pulsar nebula where a low density surroundings has delayed the action 
of the reverse shock wave. 



DISCUSSION 



Pulsars are expected to preferentially inject their wi nd power close to the equato 



rial plane and models that allow for such power injection (IKomissarov fc Lyubarskyl 12003 



Del Zanna et al.ll2004l ) have been successful in reproducing many of the properties of the torii 
observed immediately surrounding the pulsars in young wind nebulae. However, on a larger 
scale (several inner torus radii), the young nebulae such as the Crab and 3C 58 do not show 
the toroidal structure. Here we have investigated the evolutionary status of several pulsar 
nebulae, best exemplified by the newly formed Vela nebula, that do show a thick toroidal 
structure. We suggest that the structure is not present in the Crab and related objects 
because of Rayleigh- Taylor instabilities at the interfa c e between th e expanding supernova 
gas and the accelerating pulsar nebula (IChevalierlll9771 ; lHesterll2008f ) , although there remain 
uncertainties in the ins tability when the light fluid is magnetized ( IBucciantini et al.ll2004l ; 
Stone fc Gardiner! 120071 ). The instabilities give rise to mixing and diffusion of the relativis- 
ts particles. In the cases discussed here, the instability is absent either because the outer 
boundary is determined by a somewhat asymmetric outer pressure or the pulsar power has 
declined so the outer boundary is no longer accelerated. The age at this phase is ~ 10 4 yr, 
within a factor of a few. In the case of an asymmetric outer pressure, a relic pulsar nebula 
may also be present, as in the case of the Vela nebula. At later times, the pulsar velocity 
leads to a nebular trail to the pulsar. 

Multidimensional simulations are ne eded to confirm the scenario proposed here. In their 
numerical simulations of pulsar nebulae, ICamus et al.l (120091 ) find that waves and vortices in 
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the larger nebula feed back on the structure at the termination shock, which in turn generates 
more structure in the nebula; they note that it is impossible to distinguish between the cause 
and the effect. In addition, some thermal gas from the supernova is entrained in the unstable 
region, explaining the optical filaments seen in association with nonthermal filaments in the 
Crab and 3C 58. These were 2-dimensional simulations, but did not impose symmetry about 
the equatorial plane. More investigation of the instabilities and 3-dimensional simulations 
are ultimately needed (e.g.. iMizuno et al.ll201lh . 



We are grateful to Richard Dodson for permission to use the image shown in Fig. 1, and 
to the referee for helpful comments. This research was supported in part by NASA grant 
NNX07AG78G. 
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Table 1. Properties of Pulsars and their Nebulae 



PSR 


Supernova 


Distance 


E/h 5 


T~char 


Rtorus 


Rradio 


Refs. 




Remnant 


(kpc) 


(ergs s _1 ) 


(y r ) 


(PC) 


(pc) 




B0833-45 


Vela 


0.29 


6.9 x 10 36 


1.1 x 10 4 


0.030 


0.17 


1,2,3 


J2229+6114 


G106.6+2.9 


3d 3 


2.2 x 10 37 


1.05 x 10 4 


0.13d 3 


1.34 


4,3 



Re ferences - (1) 



(4) 



Kothes et al 



Dodson et al. j2003bh : (2) iDodson et al.l j2003ah : (3) iNg fc Romanil j2004h : 
J2006h 
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Fig. 1. — The image, reproduced from Fig. 10 of lDodson et al.l (j2003al ). shows the Vela nebula 
at 5 GHz, with derotated magnetic field lines and Chandra X-ray contours overlaid. The 
cross marks the pulsar position. The solid black line is the suggested contact discontinuity 
between the shocked pulsar wind and gas that has passed through the bow shock. The arrow 
indicates the direction of the mildly supersonic flow. 
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Fig. 2. — Model of a toroidal region with uniform emissivity seen in projection with the axis 
tilted 50° to the line of sight. The disk half- angle is 35° and the maximum radius is 5 times 
the inner radius. The color scale is linear. 



